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/^ • Abstract 

H : 

y^ , The transverse momentum distribution of soft hadrons and jets that accompany central 

hard-scattering production at hadron colliders is of great importance, since it has a direct 
bearing on the ability to separate new physics signals from Standard Model backgrounds. 
We compare the predictions for the gluonic /cj -distribution using two different approaches: 
resummed NLL BFKL and DGLAP evolution. We find that as long as the initial and 
final virtualities (kt) along the emission chain are not too close to each other, the NLL 
resummed BFKL results do not differ significantly from those obtained using standard 
DGLAP evolution. The saturation momentum Qs(x), calculated within the resummed 
BFKL approach, grows with 1/x even slower than in the leading-order DGLAP case. 



1 Introduction 

The high-energy behaviour of QCD amphtudes is described by the BFKL/CCFM equation [T1I2] 
which sums to all orders leading logarithmic contributions of the form (a^ Ins)". The next-to- 
leading logarithmic (NLL) corrections as{as Ins)" have been calculated in Refs. jS]. The NLL 
corrections turn out to have a large numerical effect, and therefore in order to obtain reliable 
predictions one needs to first understand the origin of the large numerical coefficients at NLL 
and then if possible to further resum the main part of the NLL contribution. 

In the past few years there have been a number of studies of the BFKL approach at NLL 
accuracy, see Refs. [H El El, which concentrate on the properties and the behaviour of the 
gluon Green function. In this paper we focus instead on the gluon transverse momentum {kt) 
distribution along the BFKL evolution chain and compare the predictions obtained using both 
the resummed NLL BFKL and the LO DGLAP evolution approaches. 

Besides the pure theoretical interest in further understanding the properties of QCD in the 
high-energy limit, there are also important phenomeno logical implications. Precise knowledge of 
the kt distribution of intermediate (i.e. accompanying) gluons is important for achieving a better 
understanding of the structure of the so-called 'underlying event' at hadron colliders such as the 
Tevatron and the LHC As has been emphasized in [7j, the fact that the transverse momentum 
of the emitted partons (and therefore of the soft hadrons and minijets) grows with energy (i.e. 
with 1/x) could cause problems due to a 'noisy' underlying event affecting the extraction of a 
clean new-physics signal (for example, inclusive production of Higgs bosons, SUSY particles, 
etc.) at the LHC. Therefore it is very important to be able calculate the predicted gluon kt 
distribution as precisely as possible. Of course the standard approach is to use parton shower 
Monte Carlos such as HERWIG and PYTHIA to estimate the kt distribution of accompanying 
partons/hadrons. However as these models are based on DGLAP evolution, they do not contain 
all the expected logarithmic contributions in the high-energy limit. It is therefore important to 
understand how and why the DGLAP and (potentially more realistic) resummed NLL BFKL 
approaches differ. Such a comparison will form a major part of our study. Also very relevant, 
is the high-energy behaviour of the saturation momentum Q^ — the transverse momentum at 
which non-linear effects (from gluon-gluon rescattering and recombination) become important 
and start to saturate the parton densities. 

The BFKL amplitude is usually described in terms of the Pomeron intercept (i.e. the 
singularity in the complex momentum j-plane) j = 1 + (^(7) which depends on the anomalous 
dimension (7) of the eigenfunction. In contrast, DGLAP evolution is usually described in terms 
of 7(ti;), i.e. the anomalous dimension is considered as a function of u;, the conjugate variable 
to X. 

In order to compare the DGLAP and BFKL predictions, in Section 2 we will consider the 
unintegrated gluon distribution /c,(x, fc*,//) — that is, the probability to find a gluon carrying 
longitudinal momentum fraction x and transverse momentum kt in a process with (hard) fac- 
torization scale /i — written in the form of a double contour integral over both 7 and tu, i.e. 
performing simultaneous Mellin transforms with respect to both the kt and x distributions. 
Depending on which of these integrations is performed first, one obtains either the standard 
DGLAP or the BFKL form. 



In Section 3 we introduce a resummation which is a modification of that proposed by the 
Firenze group [HI El- The idea is to modify the contibution of the poles at 7 = (7 = 1) 
corresponding to the normal (inverse) fc^-ordered DGLAP evolution. To be explicit, we include 
the full LO DGLAP (splitting function) contribution, allow for the appropriate 'energy' variable 
X = Q'^/s, and account for energy-momentum conservation. After these modifications have been 
made, it has been shown that the remaining part of the NLL correction does not exceed 7% of 
the original value. 

Solving the BFKL equation for the intercept 

j = l + uj{-f) = l + ^^xh) (1) 

71 

with respect to 7, one obtains the NLO contribution C„(as/u;)""'"^ to the DGLAP anomalous 
dimension 7(0;), equivalently the CriX~^a"ln'^~ (x) contribution to the (gluon-gluon) splitting 
function. It turns out that in leading-order (LO) BFKL the first two non-trivial terms vanish^, 
i.e. Ci = C2 = ^IQj. Therefore the expected BFKL corrections to DGLAP evolution are 
numerically rather weak (see |TT] for a more detail discussion). On the other hand, in terms 
of the intercept j the value of 0^(7) given by the LO BFKL equation is much larger than that 
coming from DGLAP. At first sight this looks like a contradiction. However the situation 
changes dramatically after the NLL resummation. Now the BFKL intercept becomes close to, 
and even slightly smaller than, the corresponding DGLAP quantity. Thus in any kinematic 
configuration where the initial and final virtualities (transverse momenta) are not too close to 
each other, the DGLAP and the BFKL predictions do not differ significantly (in agreement 
with the conclusions of Ref . pjl ) . 

The dependence of the kt distribution on the overall event kinematics is discussed in Section 4 
and the x-dependence of the saturation momentum Qs{x) is discussed in Section 5. Section 6 
contains our conclusions. 



2 The unintegrated gluon distribution 

We begin by considering the unintegrated (over transverse momentum) gluon distribution 
fg{x, kt,fi), in the form [T2] 

fg{x,kt,fx) = Tg{kt,ii) ' * / dzPgg{z) {x/z)g{x/z,kl) (2) 

where the survival probability T is given by 

Tg{k„f,) = exp ^- j^^ ^ ^i-*^ j^ dz' {Q{z' - A)0(1 -z'- A)z'Pgg{z') + n^P,g{z'))j . 

(3) 
Here the infrared cutoff is A = kt/{fi + h); /i is the factorisation scale; and np is the number 
of active quark flavours. When kt is close to fi we may neglect the double logarithms in the 

^Notethat Eq. ® below can be written in the form Xq = 1/7 + 2^(0) - *(1 + 7) - *(1 -7) = 1/7 + 0(7^). 



T-factor, since asln^(A;t//i) <^ 1, and thus T ~ 1. Then the unintegrated distribution is simply 
related to the derivative of the standard (integrated) gluon parton distribution function: 



fg{x,kt) 



- xg{x,fi'^) 



dln/i' 



(4) 

- ^i=kt 



Introducing the Mellin-transform variables, 7 and uj, conjugate to k^ and 1/x respectively, we 
can write fg as 

fg{x,kt)= — - / -—F{-f,uj) — (5) 
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where a = Ncag/n, Y = ln(l/a;) and r = \n{k^/Ql). F{'y,u) represents the input distribution 
at k"^ = Ql, where Qq is the starting scale for perturbative evolution assumed to be ~ 1 GeV. 
X is the resummed BFKL intercept, X = Xo + aXi, where the leading order (LL) contribution 
is 

Xo(7) = 2vl/(l)-vl>(^)_xl>(l-^) (6) 

and Xi is the NLL correction. The contour integration over u goes to the right of all singular- 
ities, while the real part of the anomalous dimension 7 is bounded by < Re7 < 1. 

In the DGLAP limit r ^ Y we close the 7 contour around the pole given by l/{uj — aX), 
leaving the inverse Mellin integral over u with 7 = 7(0;). In the BFKL case {Y ^ r) we close 
the UJ contour around the same pole l/{uj — aX) and write the result as an integral over 7, with 
UJ = uj{'y). In the latter form 

/. = r°^ |^e-^+-^i?(7) (7) 

J-ioo 2,7^1 

with uJsi^) given by the solution of the equation 

cJ-aX(7,cj) = (8) 

and 

5(7) = ^(^-"-(^)) X' = ^-^(^'") (9) 

UJ 

With this representation, the input distribution at Qq is absorbed into -8(7), which may there- 
fore be fitted to reproduce the data. 

Another possibility is to use the conventional Mellin (uj) representation for fg by closing 
the 7 contour around the nearest pole. Strictly speaking the BFKL function Xq of Eq. © 
has poles at each integer 7. The pole at 7 = corresponds to the normal twist-2 DGLAP 
contribution, the pole at 7 = 1 corresponds to inverse kf ordering {kt <^ Qq). The other 
poles at 7 = —1, —2, ...(7 = 2, 3, ...) are the higher-twist contributions (3, 4, ... gluons in the 
t-channel), corresponding to normal (inverse) kt ordering, hidden in the Reggeization of the 
BFKL gluons. However in practice we know that the higher-twist contribution at small x is 
small (see, for example, Ref. |IS|)- Therefore, to begin, we may neglect the poles at negative 
7, arguing that phenomenologically the input function F{'y,uj) is numerically small at negative 



integer values of 7. We then obtain 



|^,e^^(->+-^Z^(^) (10) 



where 



D{uj) = 7.MMH (11) 

and M[uj) = /q x'^g{x, Ql)dx is the known Melhn transform (x-moments) of the input gluon 
distribution. The first factor 7 in (fTTj) is due to the derivative in (JH), and all 7 dependent 
quantities are evaluated at 7 = 75(0;), the solution of ^. 

3 The resummed NLL BFKL intercept 

In this section we consider the resummed NLL contribution to the unintegrated gluon distri- 
bution defined in the previous section. We use the idea proposed in Ref. P, but with a small 
modification, which is rather simple and more convenient for our purpose. The crucial point is 
the fact that the major part of the 0(0;^) correction is actually contained in the 0(ti;) contribu- 
tion^. Next we note that the nearest, and most important, poles in the quantity X of Eq. (jH)) at 
7 = and 7 = 1 correspond, respectively, to the well known, normal- {kt ^ Qo) and inverted- 
{kt <^ Qo) ordered, twist-2 DGLAP contributions. Thus the leading-order characteristic func- 
tion Xq contains two parts. One is the twist-2 poles I/7 and 1/(1 — 7), a^^d the remainder is 
the higher-twist component, Xq . That is 

X, = i+-^+xr, (12) 

7 1-7 
where, from ©, we have 

XS"'^ = 2^(1) -^(1 + 7) -^(2 -7). (13) 

First, we have to modify the poles to include in the residues the full LO DGLAP splitting 
function. Because in the physical (axial) gauge both the BFKL and DGLAP LO contributions 
are given by ladder-type Feynman diagrams, the Mellin transform of the final (modified) ampli- 
tude in the u, 7 representation may be written in the same exponential form as Eqs. (I5l7ll0p . 
Thus the residue 1 in the twist-2 DGLAP pole at 7 = is replaced by the full DGLAP splitting 
function 

ijPgg{uj) = l+UjAl{uj). (14) 



2N, 



In pure gluodynamics (i.e. np = 0) Ai = —j^ + 0{iu). To account for the quark loop contribu- 
tion for Up y^ we have to replace Ai by 



^Analogously, in the d'y integral Q we cannot rule out the possibility of other singularities in the uj plane 
situated to the left of the leading pole at w = Ws. However, since in any case we cannot justify the BFKL 
approach for large |a;| ~ 1, in the present context we only keep the leading pole in ((TJ. 

^Recall that in the BFKL approach, to ^ as- 



where the first term in the brackets corresponds to the real quark two-loop contribution^, as 
indicated by the extra a^ factor. The second term is the virtual quark loop insertion in the 
gluon propagator, and corresponds to the Pqg term in the T factor of Eq. Q. Formally, the 
NLL correction is represented by the leading (u) term in Ai. Nevertheless here we keep the 
full u dependence of the LO DGLAP splitting kernels. The same procedure can be followed 
to modify the residue of the 1/(1 — 7) pole, which corresponds to the DGLAP evolution with 
inverse kt ordering. 

Another modification is necessary. In order to compare the BFKL and DGLAP predictions 
we have to write everything in terms of the DGLAP 'energy' variable x = Q'^/s. Due to the 
asymmetry in the definition of the energy scales, for the normal and inverse ordered contri- 
butions, we need to correct the contribution of the 1/(1 — 7) pole, and to replace (1 — 7) by 
(1 — 7 + cj) (see Section 3.3 and Eq. (60) of Ref. P) in the terms that are singular when 7 is 
close to 1. 

Putting everything together, we finally obtain the characteristic function 

X(7,^)=Xo + aXi, (16) 

where the LL intercept Xq is given by (jH))- The NLL contribution 

Y 7 7 1 — 'J + uj 1 — 7/ 

consists of a twist-2 part, which is shown in brackets, and a correction, — cuXq , to the higher- 
twist component (fT^ oi Xq, whose origin we now explain. Following Ref. ^1], we have multi- 
plied the higher-twist contribution (fT^ by a factor (1 — tu), which effectively accounts for the 
kinematical constraints and provides conservation of energy and momentum. At c<j = 1 the 
whole contribution vanishes. Note that the twist-2 part already satisfies this condition, since 
we use the full DGLAP splitting function. 

We have checked that, in the important region < 7 < 0.6, the approximation given by 
Eqs. fll3lll7p reproduces the known exact NLL BFKL 0{as) result to within 7% accuracy. This 
is illustrated in Fig. [T] For real 7, the continuous curve in the upper plot is the result of the 
exact calculation of Xi/Xq for np = (Eqs. (53,54) of Ref. |H], see also Refs. 0), whereas 
the dashed line comes from the approximation given by Eqs. ()13II17|) above. We need to know 
the intercept not only at the saddle point (Re 7 ~ 0.3 — 0.5, depending on the kinematics), 
but also in a region in the complex plane around this point (with |Im7| < 0.5). We therefore 
show in the lower plot of Fig. ^ the deviation, AXi, of our approximation of Xi from the exact 
result in] for various values of Im 7, normalized to the value of Xi for 7 = Re 7 with Im 7 = 0. 

Moreover, note that the approximation of Eqs. p3lll7p can be used even for rather large 
values of as, since from Eq. (jH]) we always have Ug > 0. Note that by using the above form we 
are effectively expanding in u, which here plays the role of a small parameter, rather than in 
the coupling a^ (i.e. a).^ This is the result of resumming the NLL (©(a^)) corrections to the 
leading-order BFKL/CCFM intercept. 

■^The kt integral over the second loop, in the DGLAP approximation, generates the I/7 pole. 
^With the exception of the quark loop contribution, see footnote 4. 
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Figure 1: The comparison of the approximation of Eqs. ()13II17|1 for the characteristic function 
with the exact NLL result. The upper diagram shows Xi/ Xq as a function of 7 with Im 7 = 0, 
whereas lower diagram shows the deviation, AXi, suitably normalised, of our approximation 
of Xi from the exact result for various values of Im 7. Recall that uj = dXo + a^Xi, where 
a = Sas/vr. 



The solutions of Eq. (jH)) for the leading singularity Wsi'j) are shown in Figs. El and El We 
compare the results obtained using four different approximations: 

(i) the dotted line corresponds to the well-known LO BFKL function X{'-^,uj) = Xo(7,0); 

(ii) the Double Logarithmic (DL) contribution X = I/7, which is of course the same for the 
LO DGLAP and LO BFKL cases, is shown by the dot-dashed line; 

(iii) the solution of (jHl) with the full resummed function X(7, u) given by Eqs. ()13II17|) is shown 
by the solid line; 

(iv) the dashed line shows the pure DGLAP result where we keep just the same LO DGLAP 
contribution that was included in (fT7|) . that is X = (1 + ujAi{uj))/j. 

At very small a^ = 0.01 the solid (resummed NLL BFKL) curve is rather close to the dotted 
one (LO BFKL). However already at as = 0.15 the NLL corrections significantly change the 
behaviour of 0;^ for real 7 (Figs. |2K,b,c and Figs. inK,b,c). In particular, the solid (NLL BFKL) 
curve becomes closer to the dashed (DGLAP) curve. Note that the resummed value of ujs 
depends weakly on the QCD coupling a^. The NLL BFKL solutions for a^ = 0.3(0.15) shown 
in Figs. |2b and Eb by the heavy (thin) solid lines are quite similar. 

/^From the plots, we see that in the important region. Re 7 < 0.6, the resummed NLL 
BFKL intercept is very similar to that for the DGLAP case. In fact for 7 ~ 0.3 the NLL BFKL 
value is even a little below the DGLAP intercept. However for large values of 7 > 0.6, the 
NLL BFKL curves go above the DGLAP curves, due an additional positive term from the pole 
1/(1 — 7 + C(j) which arises from the inverse fct-ordered contribution. The approximate equality 
of the NLL BFKL and DGLAP curves for 7 < 0.6 occurs because this positive contribution 
is compensated by the virtual corrections corresponding to gluon Reggeization, that is by the 
higher-twist poles at 7 = —1, —2, ... and •y = 2 + u, S + cu, ..., which give a negative contribution. 

As we go into the complex 7 plane, the resummed intercept decreases faster than that 
of DGLAP, and is closer to the original LO BFKL, since we become further away from the 
twist-2 DGLAP pole at 7 = 0. This decrease of the NLL intercept improves the convergency 
of the saddle-point integral. To demonstrate this behaviour of Us in the complex 7 plane, we 
plot in Figs. Eli and Ell the results for Im 7 = 0.5. The value of the intercept decreases with 
increasing Im 7 and, due to the higher twist contributions, we indeed see that in the BFKL 
case it decreases faster than in the DGLAP case. 

Finally, we can investigate the role of the quark loop corrections. In Fig. El we have consid- 
ered the (realistic) case of four light quarks [np = 4). In Fig. Owe show the analogous results 
for pure gluodynamics [np = 0). The inclusion of quark loop contributions evidently shifts 
the position of the singularity at Ug to the left; the absolute value of the (negative) virtual 
quark loop correction is larger than the real quark loop contribution. This demonstrates that 
the gluon spends part of the 'evolution time' (i.e. rapidity interval) in the form of a quark- 
antiquark state which corresponds to a lower intercept j < 1, that is to a negative value of 
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Figure 2: The 'Poineron' intercept, ujs-, obtained by solving Eq. (jH)), in the following four 
approximations: LO BFKL; DL; NIL BFKL and LO DGLAP; as described in the text. We 
include quark-loop corrections, with np = 4. 
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Figure 3: As for Fig. |2l but for pure gluodynamics, with np = 0. 
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4 The evolution of the gluon kt distribution 

In this section we use the formahsni developed above to study numerically the gluon kt distri- 
bution, and in particular how it changes as one moves along the evolution chain starting from a 
given input distribution at a particular value of x. We begin by introducing the Green function 
(G) corresponding to the unintegrated gluon distribution with the initial condition 5{l — k'j/k''^) 
at X = x', i.e. a fixed non-zero value of kt aX x = x' . According to the results of the previous 
section, we have 

G(r, r'; x, x') = P" ile7(-'-')+-.(^-2/'). (18) 

J-ioo 27ri 

Here r' = \n{k'^/Ql), y' = ln(l/a;') and Ug = (^^(7) is the solution of (jHJ. 

To evaluate the spread in the transverse momentum k't at some intermediate value oi y' < Y 

(that is x' > x), we write the final unintegrated distribution fg{x, kt) as the convolution 

fgix, kt) = J 0/.(^', K)G{r, r'; x, x') (19) 

and study the distribution of the integrand over k'f It is simplest to elucidate the procedure 
pictorially. The complete evolution is shown schematically in Fig. 0J We have the source at 
X = 1 and kt = Qq, and consider the evolution to the final gluon at x = exp{—Y) and some 
large value of kt, say 30 GeV. This point can be reached by different evolution trajectories, 
indicated by the shaded area in the diagram. Of course the exact form of the shaded region will 
depend on the approximation used to calculate the behaviour of the gluon along the evolution 
chain. For example, we shall see that, as expected, LO BFKL with running coupling a^, tends 
to populate the low kt domain much more than the other three approximations that we use. 

We are now in a position to study the correlated quantity of experimental interest, that is 
the distribution of intermediate gluons in a process where a high pt (or high mass) particle is 
produced, with (light-cone) momentum fraction x+ = x = exp{—Y). The scale of the process 
specifies the final kt of our gluon evolution. We have therefore written, in Eq. (jl9|) , the final 
gluon distribution as a convolution in k't of the gluon density at x' and the Green function 
fjlSp which describes the evolution from x' to the final point x. The quantity of experimental 
interest is thus the integrand, fg{x',kt)G{r,r';x,x'), which represents the k't distribution at 
an intermediate value of x'. In particular we wish to investigate how the intermediate k't 
distributions, obtained in the various approaches, differ from one another; or, to be more 
explicit, to see if the NLL BFKL result brings the LO BFKL behaviour close to the DGLAP 
result. It is convenient to plot the results as the ratio of the integrand to the final gluon density, 
fg{x', k't)G{r, r'; x, x')//p(x, kt). This corresponds to the 'normalized' distribution k'.^ dN/dk'^ of 
the intermediate gluons as a function of \nk'^, at the chosen intermediate value x' = exp(— y'). 

For the purpose of illustration we calculate the evolution of the gluon kt distribution taking 
the initial scale to be kt = Qo = 1.2 GeV and choosing the final kt = 30 GeV at F = ln(l/x) = 8. 
For the LHC energy {y/s = 14 TeV) this corresponds to a Et = 30 GeV jet with pseudorapidity 
rj ~ —2 (in the direction opposite to the parent proton), which can be observed in the central 
detector. If we set y' = 4, we can study the k't distribution of the accompanying jets with 
longitudinal momentum p^ = -EbeamC"^ ~ 55 GeV (in the parent proton direction). 

11 
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Figure 4: A schematic sketch of the 'cigar-hke' kinematic domain in x' and kt covered by the 
possible paths from the initial to the final evolution points. Note that the scales for kt and x'_^ 
are both logarithmic, as indicated by the exponential factors. Also note that the + subscript 
on x' indicates a light cone variable, which is only introduced to ensure that the fixed r] curve 
is a straight line. Figs. El and IHK^b show the gluon kt distribution at fixed values of y', whereas 
Fig. |Bt,d show the distributions at fixed rj. 
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The shape of the distribution depends on the initial and boundary conditions. Starting 
with fg = (5(1 — k'^/q^) at x = 1 (that is from the conditions corresponding to the BFKL 
Green function) we obtain almost flat distributions (see Fig. EK) in all four cases: LO BFKL, 
resummed NLL BFKL, DGLAP and DL. Note that we use the same labelling for the curves 
as in Figs. |21and|ni Using instead the input fg = 6{1 — x) at kt = Qo we obtain distributions 
that grow almost linearly with In kt, see Fig. ^jp. This is caused by the boundary condition 
fg{x, kt = Qo) = for any x < 1. Again the distributions in all four cases are close to each 
other. 

Next, in Figs. \^,d we show the results for more realistic input distributions: 

valence-like fg{x, kt = Qo) = v^(l - x)^ (Fig. Et), (20) 

and 

Pomeron-like fg{x, kt = Qo) = x'^'^l - x)^ (Fig. EH). (21) 

All the curves in Fig. El are calculated using a fixed a^ = 0.2 and np = 0. 

For the valence-like input ()20|). the initial gluons are concentrated at rather large x ~ 1 
and the results at low k't look qualitatively similar to those in Fig. EJd. For the Pomeron-like 
input (|2H) we have a bigger contribution from the small-x and low-scale kt region. Therefore 
the distributions are more close to those in Fig. EK- 

To fix the input condition at kt = qo, the usual procedure for DGLAP evolution, it is easier 
to work with the u representation, replacing the contour integration in the 7 plane (|7j) by 
the integral in the u plane (jlO|) . We have checked that to within a few per cent accuracy both 
representations give the same function fg even for the case when we keep only the leading pole in 
the denominator of the integrand (J3J. Strictly speaking this can only be true for np = 0. For a 
non-zero np there is a second pole in the 7 plane which corresponds to the second eigenfunction 
in the singlet DGLAP evolution where we have two equations — one for each of the gluon and 
quark distributions. In Section 3 we considered the gluon distribution only, and this second 
(quark) equation reveals itself as a new pole at a 7 value comparable with that (7^(0;)) of the 
leading pole. To obtain a precise solution we therefore have to keep the contribution from this 
second pole as well. However in order to simplify the discussion (and computations) in this 
section, we will consider the case oi np = only. 

One advantage of using the representation where we integrate over du is that, as in the 
DGLAP case, we can include the running of as{kp) simply by replacing the product 7s(co')r in 
the power of the exponent in (fTH)) and ([111)) by the integral 

ls{r-r') -^ f\,ico,r")dr" (22) 

J r' 

where r' is the logarithm of the initial virtuality and the anomalous dimension '-fs{uj,r") is 
calculated using the running QCD coupling as{r"). Note that in (fTUj) r' = 0, of course. The 
results for running a^ are presented in Fig. Elfor valence-like (j20j) (Figs. int',c) and Pomeron-like 
(PT|) (Figs.|nb,d) inputs. In addition to the distributions at y' = 4 (that is, at fixed longitudinal 
momentum pz ~ 55 GeV) which are shown in Figs. ^.,h, we also show, in Figs. Efcjd, the 
distributions at fixed gluon jet (pseudo) rapidity r^' = 2 in the parent proton direction. 
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Figure 5: The kfdN/dkf distributions of the intermediate gluons as a function of hikf at 
intermediate y' = 4, for fixed a^ = 0.2. We show the four distributions corresponding to using 
LO BFKL, resummed NLL BFKL, DGLAP and DL, labelled just as in Figs-ElandU The four 
plots show the distributions for different choices of the input distributions. 
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Figure 6: Plots a,b are as for Fig.|nt,d but with running a^ - also shown on these two plots (by 
thin lines) are the NLL BFKL distributions for the intermediate values y' = 1 and 6, as well as 
the y' = 4: distribution which is shown by the thick line. Plots c,d show the gluon distributions 
for a fixed value of pseudorapidity, r]' = 2, rather than fixed y'. 
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As expected, including the running a^ shifts the distributions to lower kf values where the 
coupling is larger. This is most evident for the LO BFKL curve in Fig-Eb, where the (Pomeron- 
like) input does not suppress the low kt region. After the resummation of the NLL corrections, 
the results depend less on the value of a^, and the ln/c| distribution of the intermediate gluons 
becomes broader, with the maximum shifted to a larger kt in comparison with the LO BFKL 
distribution. This is a result of the weaker dependence of the resummed intercept on the value 
of as, see Figs. |21and|ni On the other hand, we see that the form of the NLL fc^ distribution is 
very similar to that for DGLAP, since the intercepts a; (7) are almost equal in the region of the 
saddle points. 

In Figs.|nt,d we present the kt distributions of jets with a fixed pseudorapidity 77' = 2. In this 
case a small-fci jet has a smaller longitudinal momentum (that is, a larger y' = ln(l/x')). As kt 
increases we therefore move from a distribution with a larger y' to a distribution corresponding 
to a smaller y'. To illustrate this, we look at the fixed y' plots, Figs-Eti-jb, with y' = 4. On these 
plots we also show (by thin lines) the NLL BFKL distributions at two other fixed values, y' = 2 
and y' = 6. It is clear that the distributions at fixed rj' have a maximum, since the curves with a 
larger y' grow with kt, while for a smaller y' they decrease. In other words, the low kt domain is 
suppressed since at fixed 77 the gluon with a smaller kt has a smaller longitudinal momentum p^. 
In this case the low kt part of the distribution corresponds to an extreme evolution trajectory, 
where first the major part of the ln(l/x) = Y interval is used for BFKL evolution and followed 
then by an increase due to DGLAP evolution at almost constant x; see Fig. EJ Of course such 
evolution trajectories give a much smaller overall contribution than those where both 1/x and 
kt grow simultaneously leading to double logarithmic {a slia{l/x) In kt)"' contributions. 

The effect of inverse /cr ordering is more evident at fixed r]'. For the Pomeron-like input, 
the NLL distribution is broader than the DGLAP one (see Fig. Eli). However, for valence-like 
input, it is a bit narrower (see Fig. Et). The reason is as follows. The cigar-like shape of Fig. 0] 
is broader for BFKL, so we would anticipate a broader k'f. distribution. However at fixed 77', for 
low kt we enter the low x' region where the valence-like input vanishes. This has a larger effect 
on the more spreadout BFKL distribution than on that of DGLAP, with its thinner cigar which 
practically does not sample the region with k[ r^ Qq. Evidently, however, the effect is not very 
strong and the final NLL BFKL and LO DGLAP distributions do not differ significantly. 

This result is of great phenomenological importance, since it implies that for the kinematic 
regions described above, which are typical of what may be studied at the LHC, parton shower 
Monte Carlos based on LO DGLAP evolution should give a reasonable approximation to the 
predictions obtained using a full NLL resummed BFKL calculation. 

Recall that the k[ distributions have been shown in Figs. |S1 and IHl in 'normalized' form, 

dN{x',K) ^ fg{x',K)G{r,r';x,x') 
dlnk't'^ fg{x,kt) 

That is we have shown the density of intermediate gluons with momentum fraction x' = 
exp{—y') in events where the final gluon has been detected with momentum kt and fraction x. 
Therefore we have introduced the factor fg{x, kt) in the denominator. To obtain more insight, 
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we show, in Fig. [7J the actual k[ distributions 

dF{x\ k[) 



dlnk't'^ 



fg{x',K)G{r,r';x,x') (24) 



obtained in the four approximations, corresponding to Fig. IHb- This quantity, (I2H), represents 
the double inclusive cross section - the probability to observe both the final gluon (x, kt) and 
the intermediate gluon {x', k[). 

5 Determination of the saturation momentum, Qg 

Note that up to now we have been dealing with single ladder graphs and therefore have not 
taken any account of absorptive corrections (or gluon recombination). These effects are formally 
suppressed by a factor 1/R'^Q'^, where R > 1/Qo is some dimensionful parameter coming from 
the confinement domain. 

However it is well known that at very small x the parton densities grow and the non-linear 
two- or multi-ladder contributions must eventually become important. This occurs when the 
cross section for gluon-gluon recombination becomes comparable with the proton area ttR^, i.e. 
a ~ TiR^. Now at high energies the cross section behaves as a ~ x~^ / Q'^^^~"'\ Therefore to 
have a ~ tcR'^ requires 

Q^ = Ql (X a;-"/(^-^) . (25) 

This raises the question as to what value of 7 (and correspondingly ci;(7)) should one use in 

(J2SD? 

It can be shown that multi-ladder (fan) graphs becomes important for 7 > 7c, where 7c is 
given by the equation 

_ dujsilc) ^ ^s(7c) ^26) 

djc 1 - 7c ' 

This equation was first obtained in Ref. ^Hl using the 'wave front' method in the framework of 
the GLR equation. Then it was considered in Ref. ^Hl- The equation can also be justified on 
the basis of the Balitsky-Kovchegov equation (in conformal symmetric form) ^Tj d] . In the 
Appendix we show how ^I^\ can be derived using BFKL eigenf unctions. 

In order to solve Eq. (f^3j) for 7c, we return to the 7-integral representation of the amplitude 
given in Eq. ((7j). By evaluating the integral using the saddle point method, we see that for large 



r (i.e. a large interval of In fc^ DGLAP evolution) the value of the saddle point, 7sp ~ JagY/r, 
is small, but increases with the rapidity interval Y. When the position of the saddle point is 
such that 7sp < 7c, the non-linear contributions are negligible. However if 7sp > 7c, then the 
absorptive effects become large and we rapidly reach the saturation regime with fg{x, kt) ex: k^. 

The value of 7c = 0.37 was obtained in ^HI using the LO BFKL formalism. Evaluating 7c 
using the resummed NLL BFKL expressions for ujg given in Section 2, we arrive at almost the 
same value — 7c = 0.32 for ag = 0.3 — and, in fact, in the limit as —>■ the critical anomalous 
dimension 7c —^ 0.37. 

It is straightforward to obtain the x dependence of the saturation momentum Qs{x). Let 
us assume that we already reach the value oi Qs = I GeV at some Xq. With decreasing x we 
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Figure 7: The gluon distributions corresponding to Fig. ^p. Here the intermediate gluon 
distribution has not been normahzed to the probabihty to find the final gluon with transverse 
momentum kt and x = exp(— F). 
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need to increase Qs in order to keep the saddle point at 7sp = 7c. The equation which gives 
the saddle point of the integration in Eq. (jTj) is 

r = lnQH-r^ = F-^, (27) 

(27 1 — 7 

where we have used pHjl to get the last equality. We cannot determine the saddle point directly 
from Eq. (j77j) . since the value of 0;^ depends on the running of a^. Therefore we need to obtain 
a differential equation. This leads to 



d\Yi{l/x) l-7c' 

Solving this equation for 7^ leads to the x dependence of the saturation momentum Qg that is 
shown in Fig. |H1 Here we include the light quark contribution with np = 4. The labelling of the 
various curves is the same as in Fig. |21 We see that after the resummation of the NLL BFKL 
corrections the value of Qs grows much more slowly than in the LO BFKL case. Moreover, the 
NLL BFKL curve even drops below the DGLAP curve, since in the region of 7 = 0.3 — 0.4 
the value of 

^NLL BFKL ^ ^DGLAP^ (29) 

see Fig. |21 For a low value of a^, we find the expected ordering 

QhO BFKL ^ nNLL BFKL ^ ^DGLAP 

in Fig. |H|d. However at larger values of ag, we see Q^^^ bfkl i^g^omes less than Q^^^^^ , as 
shown in Fig|Ht,d. Due to the fact that the resummed value of 0;^ depends only weakly on a^, 
hiQl increases almost linearly with ln(l/a;), with 

Ql ~ a;-°■^^ (30) 

The power of x that we have obtained using NLL resummation is somewhat larger than the 
result of [Tn|. Asymptotically, at very large Qs (and as{Ql) ^ 1), the LO BFKL (dotted) and 
the resummed NLL BFKL (solid) curves become parallel. 

To determine the saturation momentum, Qg, from Fig. |S1 we need to know the value of Xq 
where, by definition, Qs = 1 GeV. If we were to take the model of Ref . [201 ; then xq would be in 
the range 0.2 x 10~^ — 0.3x 10~^, see Fig. 8 of (201. If we assume that xq = 10~^, then we see that 
the expected saturation region accessible to the experiments at HERA (Q^ ~ 2 GeV^, x ~ 10~^) 
lies outside the pure perturbative domain where one can safely neglect higher-twist effects, 
power corrections, etc. 

6 Conclusions 

The kt distribution of gluons emitted along the gluon chains that accompany a hard subprocess 
at high energy hadron colliders could give rise to two types of problem. First on the theoretical 
side, we have the possibility of infrared instability at low kt- Due to the presence of inverse 
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fci-ordering, LO BFKL evolution tends to populate the low kf domain where the QCD coupling 
as is larger. The other potential problem is more experimental. The kt diffusion, caused by the 
presence of inverse fct-ordering, may produce too many large kf gluons accompanying the hard 
subprocess, and, as a consequence, obscure the extraction of New Physics. 

To study these issues, we have used the prescription of Refs. [HI IH] for NLL BFKL resum- 
mation to calculate the resummed BFKL intercept and the k^ distribution of the intermediate 
gluons emitted during the evolution in both the BFKL and DGLAP cases. We have shown 
that after such resummation the BFKL intercept becomes much smaller than that obtained 
using the LO BFKL formalism, and rather close to that obtained using DGLAP evolution. 
Indeed, the NLL resummation tames both the potential problems. The infrared convergency 
of the distributions is similar to that for DGLAP. Secondly, for high kt, the NLL BFKL gluon 
distribution only slightly exceeds the DGLAP prediction. Thus, contrary to what might be 
expected from LO BFKL studies, we do not expect any significant increase in accompanying 
gluon (and therefore minijet and hadron) transverse momentum over and above that which 
is obtained using the parton shower Monte Carlo approach. It turns out that the presence 
of the terms with inverse ordering in kf are approximately compensated by high-twist gluon 
Reggeization. 

Moreover, it follows that the saturation momentum calculated using the resummed NLL 
BFKL framework is much lower than that predicted by LO BFKL. This implies that we can- 
not trust the numerical results obtained in the framework of the Balitsky-Kovchegov equation, 
which is based on LO BFKL. Therefore the growth of the saturation scale Qs should in prac- 
tice be much slower than that predicted, for example, in Refs. |2I1 1221 I2B|- Indeed from the 
NLL BFKL curves in Fig. |H1 we conclude that saturation effects are not expected in the pure 
perturbation region accessible to the experiments at HERA. Even at the LHC, it will not be 
easy to observe saturation in the perturbative region. Note that, from (j3Uj) . we have 

Ql ~ 1 GeV^ (a;o/x)°-^5 = 1 GeV^ (a:ov^/Q.)°-'^ (31) 

Solving this equation at the LHC energy, y/s = 14 TeV, we obtain in the centre of the pseudo- 
rapidity plateau 

Qs{v = 0) ~ 1 GeV (xov^/l GeV)°-^^ ~ 1.1 (1.6) GeV, (32) 

assuming that Xq = 10""^ (10^^). Thus, at the LHC, the only chance to observe saturation 
phenomena, with the scale not too low, Q^ ~ 10 GeV^ say, is to study processes in the 
fragmentation region, with X- ~ 0.1 and x+ ~ 10~^. 

The kt distribution of intermediate gluons emitted along the evolution chain depends strongly 
on the form of the input conditions. For a realistic 'Pomeron-like' input the distributions pre- 
dicted by the resummed NLL BFKL and DGLAP evolutions are rather close to each other. 
However, because of the contribution of configurations with inverse kt ordering the NLL BFKL 
distribution is slightly wider (i.e. has a larger dispersion). 
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Appendix: equation for the critical anomalous dimension, 

7c 

The non-linear BFKL evolution equation is^ 

dfg{x,kt 



dY 



J dki'K{k„ k[)f,{x, K) - J dk[' V ■ /|(x, k[) (33) 



where K(kt, k'^) is the BFKL kernel and V is the triple Pomeron vertex. If we expand the 
amplitude in terms of the conformal eigenfunctions E^ = e^^ (for the case of our forward 
amplitude with n = 0), 

f, = le^^Fij,Y)dj, (34) 

then we obtain the following equation for F{'~f, Y) 

^£^jX1 = uj,{^)F{^,Y) - J J V ■ F{^^,Y)F{^,,Y)d^,d^,, (35) 

where Y = ln(l/x), and where the triple pomeron vertex in the 7 representation is denoted 
by V. In addition to the extra «« (in comparison with Ug), the vertex V contains the factor 
1/R^Q^ = 1/R'^k^, which represents the small probability of recombination of two gluons of 
size ~ l/fci homogeneously distributed over a domain of transverse size ~ R. When going from 
Eq. (|HH|) to Eq. (jH3|) . we have integrated the product of the vertex V and the two eigenfunctions 
E^^ and E^^ (hidden in fg) over the transverse momentum k'^. The integral over dk[ gives a 
pole 1/(1 + 7 — 71 —72) which reflects the conservation law for the anomalous dimensions. One 
of the integrals, say dji, may be done by closing the contour on this pole, while the other (^72) 
can be performed by the saddle point method. Because of the symmetry of the expression, the 
saddle point gives 71 = 72. Thus Eq. (jH^ takes the form 

^^^^ = 0Os{k)F{k, Y) - vF\k/2, Y), (36) 



^Here we neglect the dependence on the momentum transfered through the amphtude and consider only the 
component with zero conformal spin n — (that is, we consider the forward amplitude with a flat azimuthal 
dependence). Recall that the eigenfunctions which have n ^ correspond to uJs < ^1 and therefore the non- 
linear contribution with n ^ will in any case be small, being suppressed by at least a factor as , in comparison 
with the linear one. 
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where the variables 7 and 7j are related to the 'wave vector' p^ k by 

7 = 1-A; (37) 

1+7 k . . 

71 = 72 = -^— = 1-2' ^ ' 

In Eq. ()36p. w now denotes the triple Pomeron vertex together with the measure of integration 
around the saddle point. 

To solve (jH?)|) we use the method of characteristics. In the stationary case, dF/dY = 0, we 
have 

uJs{k)F{k, Y) = vF\k/2, Y). (39) 

The solution of (jH^ may be written as F = a{k, Y)e^^^ where the factor / in the power of the 
exponent can be found by matching the function F(/c, Y) in the region of non-linear and linear 
evolution, neglecting the non-linear term in the latter. In the linear case fk = 0Js{k). Matching 
the first derivative^, we obtain Eq. (f^U]) . 

d{fk)/dk = f = ujs/k = dujs{k)/dk, (40) 

where, recall, /c = 1 — 7. 
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